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A number of methods are discussed which may serve for a treatment of electron correlations in 
solids. When the electron correlations are relatively weak like in semiconductors or a number of ionic 
crystals one may start from a self-consistent field calculation and include correlations by quantum 
chemical methods. An incremental computational scheme enables us to obtain results of high quality 
for the ground state of those systems. A number of examples demonstrates that explicitely. 

Solids with strongly correlated electrons require the use of model Hamiltonians. With their help 
one can tackle the problem of determining spectral densities for those systems. The projection 
technique is a useful tool here. In strongly correlated / electron systems electron or holes can 
crystallize with quite different physical consequences as in the case of a Wigner crystal or Verwey 
transition. Finally, different routes to heavy-fermion behavior are discussed, another hallmark of 
strongly correlated electrons. 
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I. INTRODUCTION 



The problem of electron correlations is one of the most fascinating and challenging of condensed matter physics. It 
is also of great practical importance. An example is the energy gap of a semiconductor which we certainly would like 
to be able to calculate without any adjustable parameter. Its size depends crucially on electron correlations. Therefore 
we can hope to calculate it only if we are able to treat correlations without making uncontrolled approximations. Most 
people are unaware that this goal has not been reached yet. For a detailed understanding and accurate treatment of 
correlations it is mandatory to calculate the wavefunction of the ground state and the excited states of the electronic 
system. There are other powerful methods used in electronic theory which base their advantages by avoiding just that, 
like the density functional theory (DFT) and the local density approximation (LDA) to it [1-4]. The latter has been 
extremely successful and useful from various points of view. Its conceptional simplicity is particularly noteworthy. 
However, by construction it excludes detailed insight into the correlation problem which can be gained, e.g., by 
studying the pair-distribution function g(r,rf) of the system. The latter is defined by 

flfcr') = ^_<<I>| ]T 6(S-r i )6(r-r j )\Z>. (1.1) 

Here | $ > is the ground-state wavefunction of the system, r^rj are the coordinates of the electrons and p(r) is 
their density at point r. This function yields the probability of finding an electron at point r[_ provided there is one 
at point r, relative to the one without that constraint. It is therefore suitable for describing electrons which avoid 
each other in order to reduce their mutual Coulomb repulsion. This gives rise to a correlation hole which is associated 
with each electron. An accurate description of this correlation hole, in particular at small distances is the heart of 
the correlation problem. It is for that reason that the pair-distribution function is intimately connected with the 
exchange-correlation energy E xc [p] in density- functional theory. 
The relationship is 

E xc [p] = y J d*rd 3 r>p(r) ^'J' 1 M (1-2) 
where g(r,r') is related to g(r,r') via 

l 

9(L,r') = J dXg{r,r!,X). (1.3) 
o 

The parameter A multiplies the interaction constant e 2 of the electrons, which enters g(r,r'), i.e., e 2 — > Ae 2 [5,6] when 
g(r,r') is calculated. The simplicity of the LDA or comparable approximations to density functional theory results 
from making a reasonable ansatz for g(r,r') instead of trying to calculate it from first principles. 

Given that one aims at calculating the many-electron wavefunction of a periodic solid one is faced with a number 
of quantum-chemical methods which are designed to do just that for small molecules. Particularly well known are the 
Configuration-Interaction (CI) method [7], the Coupled Electron Pair Approximation (CEPA) [8-10], and the Coupled 
Cluster (CC) method [8,9,11,12]. These methods have in common that they start out from a self-consistent field (SCF) 
Hamiltonian Hscf and its ground-state wavefunction. Correlations are included by admixing SCF excited states to 
the ground state. For solids it is important to use methods which are size consistent, i.e., which yield a correlation 
energy proportional to the system size when the systems are large. This excludes the CI method which lacks that 
property. It is a variational method and therefore has the advantage of giving upper bounds to the energy. The 
size-consistency problem can be eliminated by expressing the energy in terms of cumulants [13,14]. This ensures that 
all approximations remain size consistent. 

A SCF calculation will be a good starting point for a correlation calculation only when the correlations are not too 
strong. This will generally be the case for solids which involve s and p electrons only. When / electrons come into 
play, or in many cases also d electrons, then the electron correlations are usually so strong that a SCF wavefunction 
is not a good starting point for a computation of the exact wavefunction. In a molecular calculation the problem may 
be circumvented by starting instead from a multiconfiguration self-consistent field calculation (MC-SCF). By a proper 
choice of the active space (or of the different Slater determinants which are used) one can include the most important 
correlations already at that level and treat the remaining weaker ones by a subsequent CI or CEPA calculation. For 
a solid this is not possible because several configurations per lattice site would be required, i.e., their number would 
become infinitely large. Possibly, ways can be found in the future to solve this problem. At present, however, we 
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depend on treating strongly correlated electrons by simplified model Hamiltonians. Therefore we will subdivide here 
the discussion of correlations into one of weakly and one of strongly correlated systems. 

As regards weakly correlated systems another topic requires attention. When applying methods like CI or CEPA 
in molecular calculations, the excited states which are admixed with the SCF ground-state wavefunction arc usually 
the SCF ones, which extend generally over the whole molecule. Going over to an infinite solid this would imply the 
admixture of infinitely many excitations (configurations) to the ground-state wavefunction | $scf >■ On the other 
hand, it is evident that the correlation hole of an electron is a rather local object, which is confined to the immediate 
surroundings of it. Therefore, it should be described by local operators acting on | <J>scf > or local excitations. 
This was first consequently done in Ref. [15]. The same ideas were later applied in slightly modified versions, e.g., 
in Refs. [16,17]. Working with local operators is also useful in order to limit efficiently the number of configurations 
one is admixing to | $scf >■ This brings us to the concept of partitioning and projection which has been stressed 
by Lowdin [18]. Instead of constructing the exact ground state | ipo > from | $scf > by applying operators chosen 
from the full operator- or Liouvillc space 5i one limits their choice to a relevant subspace 3?o- The latter must contain, 
e.g., all those local operators which are necessary in order to describe the correlation hole of the electrons with a 
required degree of accuracy. Thus, instead of expanding the ground-state energy in powers of the residual interactions 
H res = H — Hscf one partitions 9? into a relevant part 3?o an d an irrelevant one 3?i and projects the wave operators 
il defined by 

I tpo > = & I $SCF > (1-4) 

onto iRo • Within that restricted space the ground-state wavefunction and energy are found by diagonalization of the 
Hamiltonian. 

The concept of projection and partitioning can also be applied to the computation of Green's functions [19]. This 
is very useful when we calculate the spectral density of a strongly correlated electron system described by a model 
Hamiltonian. The ones of Ni metal and of Cu — O planes, the most important structural element of the high T c 
superconducting cuprates serve as examples here. 

A particularly interesting phenomenon is that of electron or hole crystallization. The conditions for its occurrance 
were first determined by Wigner [20], who considered a homogeneous electron gas. A prerequisite for electron crystal- 
lization is that the electronic Coulomb repulsion dominates the kinetic energy, a condition fulfilled for a homogeneous 
electron gas only when the density is very low. In intermetallic rare-earth compounds the 4/ electrons have a small 
kinetic energy even at large densities because of the small spatial extent of the 4/ orbitals. Therefore, electron or hole 
crystallization is expected to occur in favourable cases even at high densities. There is clear evidence that YbiAsz, is 
an example of 4/ hole crystallization at low temperatures [21,22]. 

One much investigated feature of strongly correlated electron systems is a high density of low- lying excitations. 
They involve predominantly spin degrees of freedom and lead to heavy-fermion behaviour (for reviews see [23-27]). 
There is a one to one correspondence between the low-energy excitations of a heavy fermion system and those of a 
nearly free electron system, provided one renormalizes the effective mass of the latter. The characteristic low-energy 
scale of a heavy-fermion system can be of different physical origin like the Kondo effect [28] , the Zeeman effect [29] or 
quasi one-dimensional spin chains [30]. The materials CeAl 3 [31], Nd2- x Ce x Cu04 [32] and Yb^Ass [33] are examples 
for these three different physical cases. 

II. WEAKLY CORRELATED SYSTEMS 

We start out describing the formalism used for the calculation of the ground state and some of its properties of a 
weakly correlated electron system. This is followed up by a presentation of the results for a number of semiconductors 
and ionic crystals. 

II. 1. Projection techniques 

In order to ensure size consistency in the sense described above we express all size-extensive quantities in terms of 
cumulants [13,14]. 

We divide the Hamiltonian H into two parts 

H = H Q + H 1 (2.1) 

and assume that the ground state | $o > of H n is known while Hi must be treated approximately In the case of 
weakly-correlated electron systems we shall identify H n with the self-consistent field Hamiltonian H$c f and Hi with 
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the residual interactions, but we prefer keeping the discussion of the projection technique more general. We write for 
the ground-state energy 

E = (H\n) (2.2) 

where 

n = lim ( 1 + 1 - H\ \ . (2.3) 

Here the superoperator or Liouvillean L refers to H , i.e., it acts on operators A according to 

L Q A = [H ,A}_. (2.4) 

The brackets in Eq. (2.2) are defined as follows 

(A | B) = < $ I A+B | $ > c • (2.5) 

The upper script c indicates that the cumulant of that expectation value must be taken. The cumulant of a product 
of operators is defined by 



d\i" d\ n 

I — 1 



(2.6) 

A 1= ...=A„=0 



For two operators this expression reduces to 

< A 1 A 2 > c = < A 1 A 2 > - < A 1 >< A 2 > . (2.7) 



Using cumulants requires care [34]. For example, we must distinguish between the number 1 and the unit operator 

< 1 • A > c = < A > and < l op A > c = 0. (2.8) 



I op . The difference is 



^From the properties of cumulants it follows that for any operator A (but not for c numbers) the following relation 
holds 

(A\Hfl) = 0. (2.9) 

The projection or partitioning method approximates |f2) by admitting for its construction only operators from a 
reduced relevant subspace 3?o of the full operator space 3?. The operator f2 is thus projected onto 9?o- The quality of 
the approximation can be improved by successively increasing the dimension of 

Let 9?o be of dimension r and spanned by a set of operators {A v } with v = 1, r . We may choose the A v to be 
orthonormal, i.e., 

(i4„|A„) - 5^ . (2.10) 

The wave operator is then given by 

I") = I 1 + $>„A„). (2.11) 

V 

In case that all A v couple directly to Hi, i.e., when 

{A v \Hi) ^ for all v (2.12) 
one may use Eq. (2.9) to determine the coefficients r] v , i.e., 

(A^Hfl) = {A^Hi) + J2vA\\HA„) (2.13) 

1/ 

= 0. 

The solution of this set of equations can be written in the form 
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Vv = -]T L~l (A^Hi) (2.14) 

where 

L pT = (A P \HA T ). (2.15) 
The shift in the ground-state energy due to Hi is given by 

6E Q (Sl ) = J2vv (Hi\A„) (2.16) 

V 

= E *(")■ 

V 

The energy shift consists of a sum of contributions e(i/) from the operators A v which span the reduced Liouville 
space 3?o- If contains operators which do not couple directly to Hi implying that (A p \Hi) = for some A p , 

the latter enter only indirectly the calculations by modifying the coefficients r) v . In that case they have the form of a 
generalized continued fraction [13]. 

In order to give an example, we identify Hq with the self-consistent field Hamiltonian Hscf of H and Hi with 
the residual interaction part H res . We choose for the set {A u } the single (S) and double (D) excitations out of the 
SCF ground state \<&scf >• Single excitations are generated by u l = c\ where the c M , cf operators destroy and 
create electrons in the molecular spin orbitals (p^ and <fii, respectively. It is common to use Greek indices for orbitals 
occupied in \3>scf > an d Latin indices for unoccupied or virtual orbitals. Similarly, u> l ^ v describes double excitations. 
By using compound indices K and Y we can write 

A? = { ^ (2.17) 
{ ^ v ; 1 < 3 and M < v. 

The energy shift (2.16) is here equal to the correlation energy E corr and the e(v) are 

e(r) = J2 Vr(Hres\A«). (2.18) 

K 

The above choice of the operator set {A v } and the form (2.18) of the individual contributions to the correlation energy 
are identical to the coupled-electron-pair- approximation {CEP A — O) in quantum chemistry. 

It poses no problem to establish also a link between Eqs. (2.2, 2.9) and the Coupled-Cluster method which was 
originally developed in nuclear physics [12] but has been successfully applied also to quantum chemistry [8,11]. As 
has been shown in Refs. [34,35] can be brought into the form 

| ft) = |e 5 ) (2.19) 

where S is a prime operator, i.e, one which is treated as an entity when cumulants are calculated. We may 
decompose S into 

S = J>^, (2.20) 

i.e., into a basis of prime operators which have the additional property that S^l^o >^ 0. Prime operators A 
with A|$ >= need no further consideration. The i]^ are determined from (S u \He s ) = 0. Identifying the with 
the Ap (see Eq. (2.17)) we end up with the Coupled Cluster approach with single and double substitutions (CCSD). 



II. 2. Incremental method 



We can perform correlation calculations of high accuracy for the ground state of semiconductors or insulators 
applying the method of increments [36-40] . 

Let us first consider the elemental semiconductors or the III-V and II- VI compounds which all have well defined 
bonds. We express |$scf > in terms of (localized) Wannier orbitals rather than (delocalized) Bloch states. The 
Wannicr orbitals define the different bonds I. The relevant operators are identified with the operators Aj(n) and 



5 



Au(m). The Ai(n) describe one-and two-particle excitations out of bond I where n is a tcounting index. Similarly, 
the Ajj(m) describe those two-particle excitations, where the initial state of one excitation is the Wannier orbital I 
while that of the other one is the Wannier orbital J. Therefore 

= + m(n)A I (n)+ ^ m ,j(m)A I ,j{m)) (2.21) 

I ,n <IJ>m 

where the r]i(n) and rju{m) are determined according to Eq. (2.13). The method of increments provides a way of 
solving approximately that system of equations. It is based on the following idea. First, all electrons are kept frozen 
except those in bond I. After the corresponding correlation energy SE = e(I) has been determined from Eq.(2.16), 
one repeats the calculations by releasing electrons in bonds I and J. This yields a new correlation energy 5E (I, J) or 
an increment 

e(I,J) = SE (I, J) - e(T) - e(J), (2.22) 

as compared with the previous case. Thereafter, electrons in bonds /, J, K are released with the remaining ones 
kept frozen and so on. After each step a new correlation energy 5E (I,J,K) etc. is obtained with corresponding 
increments 

e(J, J, K) = SEo (I, J, K) - e (J, J) - e(J, K) (2.23) 
- e(J, K) - e(I) - e( J) - e(K) 

and so on. The total correlation energy E corr is then written in terms of the energy increments as 

E corr = Y,<I) + E < I > J )+ E <I,J,K) + - (2-24) 

/ </./> <IJK> 

where < I J >,< UK > denotes pairs and triples of bonds, respectively. The advantage of this expansion is 
that the series is usually rapidly convergent [36] even when H res is not small. It has been demonstrated that there 
are strong relations between the method of increments and Faddeev's equations, provided the latter are properly 
generalized to a hierarchy of cumulant equations [34] . These relations enable us to state the size of the error when 
the series (2.24) is terminated at some point. 

When ionic crystals are considered the increments refer to different ions instead of bonds. For example, the first 
increment to be calculated for MgO is the correlation energy attributed to an 2 ~ ion surrounded by Mg 2+ and 
2 ~ ions. An accurate modelling of the surroundings is essential, since an isolated 2 ~ is not stable. Similarly, the 
correlation-energy increment Mg — > Mg 2+ has to be determined. The two-body increments involve pairs of O — O, 
O — Mg and Mg — Mg ions and so on. 

An important point is the following. The matrix element (H res \A^ ) and (A^, \HA^) which are needed to determine 
the correlation energy (see Eq. (2.18)) can be computed from clusters. When the cluster is chosen so large that a 
matrix element does not change any more, we are sure that its value is that for the solid. This enables us to 
apply quantum chemical program packages like MOLPRO [41] to determine them and hence the different correlation 
contributions. 



II. 3. Results for semiconductors and ionic crystals 

Starting point for the computation of the ground-state wavefunction and energy is a SCF calculation. Obtaining 
accurate SCF energies for solids is still a difficult computational task. Pisani and co-workers have developed the 
successful program package CRYSTAL [42] which allows for SCF — LCAO (linear combination of atomic orbital) 
calculations on periodic solids. It expands the Bloch orbitals in terms of Gaussian-type orbitals (GTO's) and obtains 
the canonical SCF solutions. The latter can always be transformed to a localized representation by using a suitable 
localization procedure [43]. However, convergency problems are encountered when Gaussian functions with small 
exponents are used. An alternative ab-initio embedded-cluster approach has been developed within the framework 
of the LCAO method in which the localization procedure is part of the process leading to the SCF solutions [44]. 
Thereby the solid is modelled as a central cluster embedded in the field created by the remaining part of the infinite 
solid. Test calculations for LiH agree with those obtained from the CRYSTAL program, but large scale applications 
have not been performed yet. 

A third possibility for obtaining the SCF ground-state wavefunction of a solid is by means of cluster calcu- 
lations. A fragment of the lattice is used with dangling bonds saturated by hydrogen atoms. An example is 
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X2(,H^(X = C, Si,Ge, Sn). The calculated localized orbital for the central bond of the cluster can be used to 
very good approximation for every bond of the solid. 

After the SCF ground-state wavefunction has been determined the correlation corrections are incorporated either 
on a CEPA-0 (see Eqs. (2.17-18)) or CCSD (see Eqs. (2.19-20)) level, thereby applying the method of increments. 
We discuss first the results for the elemental semiconductors and the cubic III-V compounds before we turn our 
attention to the ionic crystals. 

Semiconductors 

The SCF calculations for the valence electrons are done with a basis set consisting of GTO's. Its size is 
(AsApld) I '[3s3pld] in standard notation, i.e., four s-like GTO's are contracted into three orbitals and similarly for 
p-like functions. The inner shells are described by a quasirelativistic pseudopotential [45]. In some cases a polariza- 
tion potential is added in order to account for the polarization of the inner shells by the valence electrons. For further 
details concerning the basis set we refer to the original literature [46]. The cohesive energy obtained within the SCF 
approximation is only between 55-72 % for the elemental semiconductors and between 55 - 67 % for the cubic III-V 
compounds (see Tables 1,2). Thus, even for these weakly correlated systems the correlation contribution to binding is 
appreciable. Lattice constants come out much better in SCF approximation, the errors being less than 2% (see Fig. 
3). Bulk modulili are generally overestimated by up to 20 % (see, e.g., Table 4). For a treatment of correlations the 
basis sets have to be enlarged. Instead of using one set of d functions only, we use two sets of d functions and one 
set of / functions. The basis is therefore of size (4s4p2dl/)/[3s3p2dl/] which is very respectable for diamond and 
silicon, but less for heavy elements like Ge or Sn because of low lying d states. Results for the cohesive energy of the 
group IV semiconducting and the cubic III-V compounds are listed in Tables 1 and 2, respectively. One notices that 
about 90 % of the binding energy are obtained in all cases. The missing percentage is due to basis set limitations. 
By including a second set of / functions and one set of g functions one would further improve the accuracy of the 
calculated values, in particular for the systems involving heavy atoms. 

It is well known that bond lengths or alternatively lattice constants are relatively well reproduced by a SCF 
calculation. The same holds true for the group IV semiconductors as is seen from Table 3. The inclusion of a core 
polarization potential is important in order to come much closer to the experimental values. When correlations among 
valence electrons are included in the way described above the error for the lattice constant is less than 0.7 % in all 
cases. 

The bulk modulus B is defined by 



where Eo(a) is the ground-state energy as function of the lattice parameter a. As mentioned before, self-consistent 
field calculations overestimate generally B. The core-polarization potential as well as the correlation energy lead to 
a decrease of B so that the final results are close to the experimental ones (see Table 4). This is easy to understand: 
when correlations are taken into account they keep electrons apart. Therefore an increase in pressure results in a 
smaller energy increase than obtained by a SCF treatment, where electrons come too close to each other at the expense 
of Coulomb repulsion energy. 

We are particularly interested in trends concerning the interatomic correlation energy. The latter is obtained 
when a minimal basis set is used, i.e., one containing one set of s and p functions only. Correlations, which require 
for their description a larger than minimal basis set are termed intra-atomic [14]. Interatomic correlations reduce 
charge fluctuations between neighboring atoms. The simplest way of obtaining a good estimate of the one-bond and 
two-bond contributions is by making a bond-orbital approximation (BOA). In this approximation the ground state 
is written as a product of independent bonds containing two electrons each. In that case the interatomic correlation 
energy depends on the following quantities: the interaction matrix element Vq describing the difference of the Coulomb 
repulsion when two electrons are situated on the same atom and when they are placed on different atoms of the same 
bond; the electronic hopping matrix element to between atoms of the same bond; the van der Waals interaction, in 
simplest approximation between neighboring bonds, V\; the bond polarity a in case of heteropolar bonds [52]. 

For the elemental semiconductors the bare hopping matrix element to scales with the tight-binding band width of 
the system. For diamond the value is t$ = lO.leV as obtained from band-structure calculations. We can determine 
the one for Si, Ge and a — Sn according to Froyen and Harrison [53] from the width of the LCAO bands. Once to is 
known the parameter Vq and V\ follow from the one-bond and two-bonds increments 




(2.25) 




) 



(2.26) 
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(I and J label nearest-neighbor bonds) when the latter are computed with a minimal basis set. A list of the 
calculated parameters is given in Table 5. One can use them in order to derive a simple scaling relation between the 
size of to, Vo, V\ and the bond length d. One finds that the following relations describe the data fairly well [54] 

t°/t% = (d c /d) ie (2.27) 
V /V C = (dc/d) 1 - 7 
Vt/Vf = (d c /d) 1A 

where the script C refers to diamond. Going over to the polar III-IV compounds one expects a relationship 

E%%r(a) = (l-a 2 )^™* e /(0) (2.28) 

between the interatomic correlation energy in the presence of the polarity E™l% r {a) and the one calculated from 
t ,V~o,Vi without taking the polarity into account, i.e., £^™*^ r (0). In fact, the BOA predicts that v = 5/2. But, 
when one calculates S™*^ r (0) from Eqs. (2.28) using the experimental values of d and compares it with the ab initio 
results using minimal basis sets one finds a value of v = 4.0 ± 0.25. The difference between the two approaches is the 
following. In the BOA the excited electrons remain in their original bonds while in the ab inito approach excitations 
into other bonds are taken into account too. A simple estimate of the zniraatomic correlation energy is more difficult. 
Using the values of the free atoms is a too rough method. One way is to apply an atoms in molecules approach [55] 
originally suggested by Moffitt. It takes into account that the ratio of s and p electrons differs from that of a free 
atom when a solid is formed. Also, the electron number fluctuates when the atom is part of the solid and this is also 
accounted for in that approach. For more details we refer to the literature [14]. 



Ionic crystals 

The calculations for the ground state of semiconductors described above can be extended to ionic crystals. Also 
here, correlations contribute approximately 1/3 of the cohesive energy. We demonstrate this by studying MgO, CaO 
and NiO. Thereby we employ again the method of increments, but in distinction to the semiconductors considered 
before, the increments refer here to ions instead of bonds. SCF calculations for the binding energy, lattice constant 
and bulk modulus were done by applying the program package CRYSTAL. Results are found in Tables 6 and 7. 

Correlation calculations are performed on clusters of one, two and three ions which are embedded in their proper 
surroundings. The X 2+ (X = Mg, Ca, Ni) ions are described by small-core pseudopotcntials and extended basis sets 
or a large-core pseudopotential together with a small basis set and a core-polarization potential. Test calculations 
for single ions show that the potentials describe very well the results of CCSD calculations with all-electron basis 
sets. The latter are in good agreement with the experimental values for the first and second ionization potentials 
X — > X + ,X + — ► X 2+ , in particular when they arc supplemented by a pcrturbational treatment of the triplet 
excitations (CCSD(T)). For NiO we applied quasi degenerate variational perturbation theory (see Ref. [60]) for 
technical reasons: the MOLPRO program package does not allow for CCSD calculations of NiO with low spin. 

The one-body increments are attributed to an 2 ~ and a X 2+ ion. Consider first the 2 ~ ion. It must be 
embedded in a proper surroundings, since a free 2 ~ ion does not exist. Stabilization is achieved describing the 
six nearest neighbors by means of a pseudopotential. The crystal environment of the resulting seven-atom cluster is 
represented by 336 ions with point charges ± 2 surrounding this cluster in form of a 7x7x7 cube. Employing a basis 
set [5s4p3d2f] one obtains the increment O 2 ~ shown in Table 8. It is noticed that it is nearly independent of 
the system considered. This is different for the increment X — > X 2 ~ which increases from Mg to Ca to Ni (see Table 
8). This can be understood by realizing that for Ca, excitations into the low lying d orbitals are important while 
for Ni the correlation energy is further enhanced by near degeneracies of d orbitals. A more detailed picture of the 
one-body increments for NiO is shown in Fig. 1. For MgO and CaO the results are qualitatively similar. 

Of particular interest are two-body increments. Hereby we have to distinguish between O — O, X — O and X — X 
increments. When the distance between the two sites becomes large, then in all three cases the increments describe 
van der Waals interactions of the ions. We show in Fig. 2 the — increments for MgO for different distances. They 
do not contain their respective weight factors. Table 9 lists the different two-body increments for the three oxides. 

Van der Waals interactions can be estimated with the help of London's formula. It states that the lowering of the 
energy of two sites A and B due to van der Waals interactions is given by [62] 

"°-f> < 2 ' 29 > 

where Ia(b) are the ionization energies, (xa(b) are the polarizability volumes, R is the distance between sites A and 
B and n is a constant of order unity. It relates a proper atomic mean-excitation energy to the ionization energy. In 
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Fig. 3 we have plotted the calculated values for the van der Waals interaction between positive and negative ions 
(circles) and between the negatively charged ions (crosses) when quantum-chemical methods are applied (here CCSD) 
and when London's formula is used. The data can be roughly approximated by a straight line of slope 4. This suggests 
that Eq. (2.29) should be multiplied by a correction factor of order 4. 

Three-body increments are very small. For the three oxides the sum of them is of order 0.03 eV (MgO) 0.05 
eV (CuO) and 0.11 cV (NiO). In Fig. 4 a summary of the different increments to the cohesive energy of MgO is 
given. Table 6 summarizes the results for the cohesive energy for all three oxides. Table 7 does the same for the 
calculated lattice constants. It is found that correlations have two effects on the latter. The van der Waals interactions 
lead to a reduction of the lattice constant since the corresponding energy gain increases with decreasing distances. On 
the other hand, the correlation energy of 2 ~ increases with increasing lattice spacing since the excitation energies 
are lowered in that case. Therefore, an increase of the lattice constant is favoured by that effect. This is opposite to 
what one would expect from the LDA because the density decreases as the lattice constant increases. We find that 
the calculated lattice constants deviate by less than 1 % from the experimental values. 

Ground-state calculations based on quantum chemical methods are certainly more costly than those based on 
approximations to density-functional theory. However, they allow for systematic improvements by using larger basis 
sets. For the systems discussed here the quality of the results is certainly comparable with that of the generalized 
gradient corrected LDA. 

In addition to the ground-state calculations one would like to be able to calculate the energy band for a system like 
NiO. For work in this direction see, e.g., Refs. [64, 65] and further references cited there. 



In solids with d- or /- electrons the electron correlations are generally strong. In that case a SCF calculation is not 
a good starting point. This holds particularly true for 4/ or 5/ systems, but also a treatment of the electrons in the 
Cu — O planes of the high-T c superconducting cuprates should include from the outset the strong on-site Coulomb 
interactions of holes in the Cu id shell. Calculations of this kind can be presently done only with simplified model 
Hamiltonians. They all have in common that the basis set is a minimal one, or even less, like in the case of a one-band 
Hubbard model, which has been proposed for a study of d electrons. From quantum-chemical calculations on small 
systems it is known, that a minimal basis set gives unsatisfactory results in almost all cases. Nevertheless, the hope has 
been that an understanding of the model systems will yield generic properties of strongly correlated electron systems 
which should prevail when eventually a more realistic description should become possible one day. This applies, in 
particular, to the low-lying excitations in those systems. It should be pointed out that the local-density approximation 
to the density-functional theory does not have basis set problems. Since the contributions to the energy are calculated 
from an exchange-correlation expression for a homogeneous electron gas, characteristic features of strongly correlated 
electrons as obtained, e.g., from a Hubbard Hamiltonian are not present in that approach. This changes partially 
when a LDA+U scheme is used [66], where an on-site repulsion term U is added by hand to the systems energy. 
First, we describe the Hamilton operators most commonly used. Afterwards, we use them in order to calculate the 
spectral density for Ni with special emphasis put on the reproduction of the observed satellite peak in photoelectron 
spectroscopy [67], and of doped Cu — O planes which are part of the cuprates. 



The most frequently used model Hamiltonian for strongly correlated electron systems, is the one first suggested 
independently by Gutzwiller [68], Hubbard [69] and Kanamori[70] and commonly referred to as Hubbard Hamiltonian. 
It is of the form 



The first term denotes the kinetic energy with < ij > referring to a pair of nearest-neighbor sites. The second term 
refers to the on-site Coulomb repulsion with rii a = af^a^. For strongly correlated systems the ratio t/U <C 1. In that 
limit the Hamiltonian can be transformed into one acting on a reduced Hilbcrt space from which all configurations 
containing doubly occupied sites i are excluded. This effective Hamiltonian is of the form [71] 



III. STRONGLY CORRELATED ELECTRON SYSTEMS 



III.l. Model Hamiltonians 




(3.1) 



<ij>a 




(3.2) 



<ij>a 



<ij> 
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and is called t — J Hamiltonian. The af a , di a operators act on the reduced Hilbert space only and are given in terms 
of the af a , aw operators as 

& = o+(l-n i _ ff ) (3.3) 

The spin operators are = (1/2) XL/3 ^ta^-ap^ip where a refers to the Pauli matrices. Furthermore, hi = J2 a ^tr®™ 
and the coupling constant is J = At 2 /U. For transition metals a five-band Hamiltonian describing the rather strongly 
correlated d electron is often used. It is of the form 

H = H Q +Y / Hi(£) (3.4) 

e 

H = '^2 e v {k)n va {k) 

vak 

= \ E J2 V ^nat^)a+ a ,(£)a na/ (e)a ja (£) 

ijmn <j<j' 

where t is a site index and i, j, m, n label different d orbitals. The e v (k) [y = 1, 5) describe the energy dispersions 
of the canonical d bands. Furthermore, n va (k) = c+ a {k)c v<T (lt) The creation operators c+ CT (fc) of the Bloch eigenstates 
are expressed in terms of the af a (£) as 

4M) = -4= E (3-5) 

(N is the number of sites). The interaction matrix elements are of the form 

Vijmn — UimSijSmn + J{j {Si n Sj m + &ivn&jn) (3-6) 
Uitn = C/ ~t~ 2t/ 2(/^ m 



where f and J are average Coulomb- and exchange-interaction constants. For a cubic system the matrix can be 
expressed in terms of the average exchange constant and a single anisotropy parameter A J only. The explicit form 
of the matrix is found, e.g., in Ref. [72]. 

For a description of electrons in the Cu — O planes, the crucial structural element of the high-T c cuprates, a three- 
band Hubbard model is usually used. It takes into account only the Cu 3d x 2_ y 2 and the O 2p x (y) orbitals (see Fig. 5). 
Their orbital energies are td and e p . Two holes on a Cu or O site interact with an on-site Coulomb matrix element Ud 
and Up, respectively. The hopping matrix clement between a Cu orbital and an O orbital is denoted by t p d- Values 
for the different parameters can be obtained from constrained LDA calculations [73]. There is general consent that 
the following values are reasonable: Ud — 10.5eV, U p — 4.0eV, t p d — 1.3eV, t pp — 0.65eV, e p — e d = 3.6eV (in hole 
representation) . 

Written in a basis of O orbitals which is diagonal with respect to O — O hopping, the three-band Hubbard Hamil- 
tonian reads 

# = E e ™(k) P + kaPm k* + u p j2 n Pi( j ) n Pi( j ) + ( 3 - 7 ) 

mka J 

+ e d^dl a d ka + U d ^n d] (I)n di {I) + 

ka I 

+ 2t pd ^ {$™kV m ka d hP + (t>* mk dt a P m ka)- 
mka 

Here / and J are indices for the Cu and O sites, respectively. Furthermore, n pa (J) = Pj a p.ja and n d a(I) = dj a dia- 
The e m (k) are given by 

e m (k) = e P ± 2t pp [cosk(r 1 + r 2 ) — cosk^ — r 2 )] (3.8) 
(m = l,2) 

with vectors r 1 , r 2 pointing from a Cu site to the two O sites of the unit cell. The phase factors 4> m k are 
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4> m k = —/= [sin kr 1 ± sin kr ] . (3.9) 
v2 

The model Hamiltonians (3.4) and (3.7) are used in the following in order to calculate the spectral density of Ni and 
of the Cu — O planes, respectively. In both cases the Coulomb repulsion exceeds the hopping-matrix elements and 
therefore a simple perturbation expansion in terms of the interactions is not sufficient. Instead, the large interactions 
must be taken into account more accurately. This goal is achieved by applying projection techniques. 



III. 2. Spectral densities 

Spectral densities can be measured either by photoelectron emission (PES) or inverse photoemission spectroscopy 
(IPES). They are calculated from the positions and intensities of the poles of the single-particle Green's function. 
One way of determining those poles is treating numerically small clusters, e.g., by means of the Lanczos method 
[74]. Thereby periodic boundary conditions are assumed. Here we shall proceed differently. Using projection and 
partitioning techniques we are able to do the calculations analytically to a large extent with some numerical work 
remaining. The latter is by far less than the one in the numerical approaches mentioned before. 

We start out describing the projection technique as applied to the calculation of spectral densities. Consider a set 
of wave- number dependent operators A m (k) in Heiscnbcrg representation, i.e., with a time dependence given by 

A m (k,t) = e «(*-^) A m {k)e- lt{H -^ N \ (3.10) 

Here H denotes the Hamiltonian of the system while /x is the chemical potential and N is the operator of the total 
electron number. Using a grand-canonical ensemble and working with a fixed value of u has the advantage that the 
calculations can be readily extended to finite temperatures. 

We define a matrix of retarded Green's functions G mn {k, t) for the set of operators A m (k) through 

G mn (k,t) = -i6{t) < VoP+&i),A„(fc)] + |Vo > (3.11) 

where -00 is the ground state of the system described by H while 9(t) is the step function which equals 1 for t > 
and for t < 0. Introducing the bilinear form 

(A\B)+=<M[A + ,B\ + \^> (3.12) 

we can write the Laplace transform of Eq. (3.11) in the form 

G mn (k,z) = (A m {k)\-^—A n {k))+. (3.13) 
Z — L 

The Liouvillean L acts on operators A according to 

LA = [H,A}_. (3.14) 
The spectral functions A mn (k, uj) belonging to the set of operators {A m (k)} are defined according to 

A mn (k, u) = - — lim Im{G mn (k, u + if])}. (3.15) 

Equation (3.13) is in a proper form for applying projection and partioning techniques. The idea of partitioning 
was put forward by Lowdin [18] and applied in particular by Pickup and Goscinski [75] and Linderberg and Ohrn 
[76]. Projection techniques were also introduced by Mori and Zwanzig for the description of dynamical correlation 
functions. Here we will use an extension of the technique to static quantities in a form which is size consistent [13,14]. 

We proceed as follows: We add (or remove) an electron in a Bloch state of the system for which we want to compute 
the spectral density. One of the A m (k) is identified with the corresponding creation (annihilation) operator. The 
strong correlations are taken into account by a proper choice of local operators to which the added electron (hole) 
couples strongly. They modify the surroundings of the electron added to the probe (IPES) or of the hole left behind 
by an ejected electron (PES). We identify the most important of these processes from case to case and, after a Fourier 
transformation, include the corresponding operators in the set {A m (k)}. Examples are given later when the theory 
is applied to Ni and the Cu — O planes. All other microscopic processes are discarded. The operator space 3? is 
partitioned this way into a relevant part 3? spanned by the {A m (k)} and an orthogonal, irrelevant part 9?i which 
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is neglected. The desired spectral density is then obtained from one of the diagonal elements of the Green's function 
matrix G mn (k,u> + irj) (see Eq. (3.15)). We rewrite Eq. (3.13) in matrix notation as 

G(k,z) = ^©[^(fc)-^(fc)]- 1 ^) (3.16) 

with the susceptibility and frequency matrix defined by 

Xmn(k) = (A m (k)\A n (k))+ (3.17) 
u) m n(k) = (A m (k)\LA n (k))+ , 

respectively. There is no memory matrix appearing in Eq. (3.15) because we neglect 

The matrix elements (3.17) can be evaluated in two different ways. One consists in using a general relationship 
between spectral functions and static expectation values. For the present case it can be written in the form 

<MAUk)A n (k)\%> = J <kjA mn {k,Lj)f(Lo). (3.18) 

— oo 



Here f{ui) is the Fermi function which can be replaced by a step function. With the help of this relation the static 
expectation values (3.17) can be determined self-consistently. The second way of evaluating them is by applying the 
projection technique also to these static quantities. Thereby one uses the following relations 

< ^o\[A+(k),A n (k)} + \^ > = m<(k),A n (k)] + Q) (3.19) 

< MA+(fc),£A„(fc)] + |^ > = (mA+(k),(LA n (kyy] + n) 

which are derived, e.g., in Ref. [14]. The round brackets and il are defined by Eqs. (2.5) and (2.3), respectively. 
The notation {AB)' implies treating the product AB as an entity when the cumulant is calculated. This completes 
the description of the projection technique as applied to Green's function calculations. When the theory is applied, 
a proper choice of the operator set {A m (k)} and of ft is crucial. 



III. 3. Application to 3d transition metals 



We apply the above theory in order to calculate the direct and inverse photoemission spectra of 3c? transition metals 
with special reference to Ni. It has been known for a long time that electronic correlations influence considerably the 
excitation spectra of these materials. A well studied case is that of Ni, where angular-resolved photoemission data 
reveal a 25 % reduction in bandwidth as compared with LDA bandstructure calculations and also the appearance 
of a satellite structure 6 eV below the Fermi energy [77]. A number of different investigations have dealt with these 
experimental findings. Starting with the work of Kanamori [70] , different methods have been applied by Pcnn, Liebsch, 
Igarashi, Roth, Hertz and Edwards and others [78-82]. We will show that the projection technique is here a valuable 
tool leading to satisfactory results for Ni [83]. 

We start from the Hamiltonian (3.4) for the 3d electron and identify the e u {k) with the canonical d bands obtained, 
e.g., from a LDA calculation. The spectral density is obtained from the retarded Green's functions 

GW(fc,t) = -*6(t) < M<w(£,t),c+ (fc,0)]+|Vo > • (3.20) 

We project again onto the relevant part % of the Liouville space spanned by a set of operators {A m (k)}, one of 
which is c+ CT (fc). The associated Green's function matrix is again of the form of Eq. (3.13) which is rewritten as in 
Eq. (3.16). The susceptibility and frequency matrix are of the form of Eq. (3.19), i.e., 

Xmn(k) = (n\[A+(k) 7 A n (k)] + n) (3.21) 

u mn {k) = (n\[A+(k),(LA m (k)y}+n). 

The ground state | $ > on to which ft is acting is here the nonmagnetic SCF ground state of H, i.e., 

i^o > = n c -©i°>- ( 3 - 22 ) 

V(7,\k\<kF 
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The strong correlations are taken intoaccount by including in the set of {A m (k)} the Fourier transforms of a number 
of local, on-site operators. They are 

A^U) = / 2 <^^ (3 23) 

13 \ a+(£)S nj (£), 1^3 

The notation 8n ia (i) = n ia (l)— < n ia (i) > and s^t) = (l/2)T, a /3af a (£)a a/i aip(£) has been used here. The selected 
A m (k) consist therefore of A^\k) = c^(k) and 

4 } ©-^E4 r) w e - f ( 3 - 24 ) 

r =1,2,3. 

For a given value of k the total number of relevant operators is 66. This requires the diagonalization of a 66 x 
66 matrix for each k point. In order to evaluate the matrix elements we must first specify the wave operator ft. In 
accordance with Eq. (2.11) we make the ansatz 

|0) = |l + ^ % ^(£)) (3.25) 



where the local operators 60ij(£) are given by 



26n^{£)dn il (£), i = j 

5 ni (£)8n 3 (£), i^j (3.26) 

The operators s^tjs^t) generate Hund's rule correlations in the ground state of a transition metal, while the other two 
types of operators reduce fluctuations of charges from their average values, thus rendering the Coulomb repulsions less 
effective. More details are found, e.g., in Rcf. [14]. The resulting spectrum for Ni is shown in Fig. 6. All parameter 
values arc in units of the SCF bandwidth W. They are obtained by fitting the measured multiplet structure of 
transition-metal ions embedded in simple metals [84] . Shown is the SCF density of states and the modifications which 
are obtained when the correlations are included. In order to bring out more clearly the changes caused by the different 
interactions we show in Fig. 6 c the special case where the exchange constant J as well as the anisotropic exchange 
parameter AJ have been set equal to zero. When J = A J = one obtains only one quasiparticle and one satellite 
peak for each k point and given band index v. The satellite peak in the density of states reflects the e g — ti g splitting 
caused by the ligand field. When J is included, but A J = new spectral density appears near -2.1 and -0.7. This 
reflects the atomic d 2 (hole) multiplet which splits into a X S state, two degenerate singlets 1 G and 1 D and two 
degenerate triplets 3 P and 3 F. The splitting energy between 1 S and 1 G is 5 J while the one between 1 G and 3 F 
is 2J. The spacings between the structures at -2.1, -1.1 and -0.7 are of comparable size. Inclusion of AJ ^ splits 
the main satellite at -1.1 into smaller peaks. We also show in Fig. 6c the modifications in the spectrum which arise 
when we set = 1, i.e., when we neglect the correlations in \tpo > and replace it by the SCF ground state |$o >• 
In conclusion, we may state that the correct position of the satellite peak as well as a band narrowing by 15 % are 
obtained when experimentally determined parameters for the interactions are used. The band narrowing is less than 
the observed one which is approximately 25 %. 

The spectral density of Ni can also be calculated by using Faddeev's equations [85]. The results resemble the ones 
presented here but they are not identical [86]. Faddeev's equations also contain three-particle correlations which have 
not been included here. A serious deficiency of the model is the neglect of the 4s band. This might explain some 
discrepancies which arise when the present theory is applied to other 3d transition metals. For example, we also obtain 
multiplet structured satellites in the spectra of Go and Fe when we use the Coulomb parameters as determined from 
optical experiments. However, such structures have not been observed in photoemission experiments. Further details 
can be found in Ref. [83] . 

The long-term goal is certainly an extension of the theory to larger basis sets. One would also like to take system- 
atically into account correlations between sites, an extension not yet tried for the 3d transition metals. In principle, 
the projection technique is a very suitable tool for modern computing. 
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III. 4. Spectral functions of Cu — O planes 



The Hamiltonian for the electrons in the Cu — O planes of cuprates has been discussed before (see Eq. (3.7-9)). 
The main task for obtaining the spectral densities is the proper choice of the set of relevant operators {A m (k)}. 

The strong electron correlations, which must be accounted for, determine the selection of these operators. They 
must include, first of all, the hole operators 

A p (m,k) = p+ kv A d {k) = d+ (3.27) 

but also the Fourier transforms of a number of local operators, which generate the correlation hole surrounding 
the added particle or hole. To those belong the Fourier transforms P^fcT'^T of the local operators p^ = pf^n P i(I) 
and dj a = dj^ndi(J), respectively. They ensure reduced weights of configurations with doubly occupied Cu and O 
orbitals, a consequence of the Coulomb repulsions Ud and U p . Additional microscopic processes to be included arc 
defined by the operators [87] 



Mk) = ^e-^pt lS i 



Mh) = -!= £ e-^p+pj.dn. (3.28) 

The N different unit cells have been labeled by / and their lattice vectors by Rj. The operator Sf~ = df^dii 
describes a spin flip of a hole in the 3d x 2_ y 2 orbital of site /. The operator p^ refers to the combination 

p+ = - p+ - p+ - p+) (3.29) 

of the four O orbitals surrounding Cu site I. The operator Af(k) ensures inclusion of processes, which result in the 
formation of a singlet state between a hole on a Cu site and another one in a nearest-neighbor O orbital (Zhang-Rice 
singlet [88]). The operator A a (k) does the same for the triplet state while A c (k) takes charge transfer processes 
in the vicinity of the added particle (hole) into account. With this choice of a total of nine relevant operators a 9 
x 9 matrix G mn (k, u) must be diagonalized for each k point. Thereby the self-consistency condition (3.18) is used 
throughout. The resulting densities of states for the case of half-filling (corresponding to LC12CUO4) and of hole 
doping (e.g., La2- x Sr x CuOiC) agree very well with those obtained from numerical diagonalization of a cluster of four 
units of CuOi [89,90] (see Figs. 7,8). This proves the usefulness of the method. As seen from those figures the 
half-filled system is insulating. The structure around 2.5 t p d results from the singlet state, while the one at Bt pd 
represents the upper Hubbard band. When the system is doped (see Fig. 8 b) spectral density is moved from the 
upper Hubbard band to the region close to the Fermi energy. The system is then metallic. In electron-doped systems 
like Nd2- x Ce x CuOi the Fermi level is in the upper Hubbard band. In that case spectral density is shifted from the 
singlet states to the upper Hubbard band which can accommodate an increasing number of electrons as the doping 
increases. When the hole concentration approaches zero, the upper band has a total spectral weight of 2 as compared 
with the half- filled case where this weight is 1. For more details we refer to Ref. [87]. 



IV. ELECTRON CRYSTALLIZATION 



In metals the kinetic energy of the electrons is usually more important than their mutual Coulomb repulsion. This 
is due to Pauli's principle which results in Fermi energies of the order of a few eV. However, this changes when 
the electron concentration is very low. Let r denote the mean radius of the volume an electron has available, i.e., 
define r Q via (4Tr/3)r^ = p^ 1 where p is the electron density. The average kinetic energy Sekin of an electron due 
to the uncertainty relation is Sekin = (Ap) 2 /2m ~ 1/(2t7?Tq). But the average Coulomb repulsion is Se po t ~ e 2 /ro 
and therefore larger than Sekin in the limit of large tq or low densities. By considering a homogeneous electron gas 
with the positive charge uniformly spread over the system (jellium model) Wigner [20] discussed the form of the 
ground state in the low-density limit and concluded that electrons crystallize in form of a lattice in order to keep the 
Coulomb repulsion as low as possible. The kinetic energy reduces to the zero-point motion of the electrons around 
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their equilibrium position. The change from a homogeneous to an inhomogeneous electron charge distribution takes 
place at a value of approximately ro/as ~ 40 — 100 where as is the Bohr radius. At those densities the electrons 
are sufficiently far apart that the exchange plays only a minor role because of its exponential decrease with distance. 
The excitations are given by the vibrations of the electrons around their lattice positions and the low-temperature 
specific heat C ~ T 3 like for a phononic system. Possible realizations of Wigner crystal are found in semiconducting 
inversion layers [91, 92], or rare-earth pnictides with low carrier concentrations [93]. For a comprehensive review see 
Ref. [94] . Often is is more appropriate to consider electrons on a lattice. In order to explain the dramatic temperature 
dependence of the resistivity of magnetite {Fe^Oi), Verwey [95] developed a model for charge ordering in that rather 
complex spinel structure. According to his model the 2 • Fe 3+ + 1 • Fe 2+ ions per formula unit are distributed at 
low temperatures as follows. One Fe 3+ is used to form a stable cubic sublattice. The remaining Fe 3+ and the Fe 2+ 
ions form a structure on which they alternate. This is achieved by assuming that the Fe 3+ — Fe 3+ nearest-neighbor 
interaction energy is higher than the Fe 2+ — Fe 3+ one, so that Fe 3+ ions prefer Fe 2+ ions on nearest neighbor 
sites. The kinetic energy of the electrons is discarded in Verwey's theory. According to this model a first-order phase 
transition of order-disorder type takes place at high temperatures. The ground state suggested by Verwey can be 
considered as a form of electron crystallization which is quite distinct from the one considered by Wigner. Its explicit 
form has been questioned though, by neutron-scattering experiments. 

A third form of electron crystallization is due to Mott [96,97] and Hubbard [98,69], and is also based on electrons 
positioned on a lattice. Here it is the on-site Coulomb repulsion of electrons which may lead to crystallization, 
provided the system is at half-filling (i.e., with one electron per site), and the Coulomb repulsion is sufficiently large 
as compared with the hopping matrix element between sites. Mott realized that this is always the case if a chain of 
H atoms is considered and the latter are pulled sufficiently apart. Hubbard discussed the metal- insulator transition 
associated with electron-crystallization by suggesting various approximations for the computation of the excitation 
spectrum of the Hamiltonians (3.1) at half filling. 

Recently it has become clear that still another modification of electron crystallization is realized in Yb^Asg, 
[30,22,99,100]. There have been previous observations on charge ordering on Yb^Ass and SmiBis, without in- 
terpretations offered [21, 101]. Like in the case of the Wigner lattice it is here the long-range Coulomb interaction 
which results in what is actually a crystallization of holes. But in distinction to the Wigner case crystallization takes 
place at high densities. This is so since 4/ holes are involved here which have a very small kinetic energy due to the 
small hybridization of the well localized 4/ orbitals. 

Yb^Ass is of the anti-T/i 3 P4 structure. The Yb ions are positioned on four families of interpenetrating chains which 
point along the four diagonals of a cube. It is important that the distance of neighboring Yb ions within a chain is 
3.80 A and therefore larger than of Yb ions between different chains which is 3.40 A. The three nearest neighbors of 
a Yb ion therefore belong to the other three families of chains. The structure is shown in Fig. 9. 

We write Ybl + Yb 3+ As 3 in order to demonstrate that the system has one 4/ hole (Yb 3+ ==>• 4/ 13 ) per formula unit. 
At high temperatures the 4/ holes are delocalized and the system is metallic. Measurements of the Hall coefficient 
show that the carrier concentration is indeed approximately 1/4 per Yb ion. At T$ = 300 K the system undergoes 
a weak first-order phase transition below which the Yb 3+ ions accumulate on one family of chains, e.g., those along 
the [111] direction. Since the Yb 3+ ions are smaller than the Yb 2+ ions, the phase transition is accompanied by a 
volume conserving trigonal distortion. Thereby chains in the [111] direction are shortened while those in the other 
three directions parallel to the diagonals of a cube are elongated. 

The driving mechanism of the phase transition is the Coulomb repulsion of holes. The energy is minimized when 
the holes move into one family of chains because of the large distances between ions in a chain. Disregarding first the 
structural changes associated with the phase transition we are faced with a Hamiltonian for the holes of the form 



The operators af^, dj a are the same as in Eq. (3.3). The first term describes the kinetic energy of the 4/ holes 
and the prime indicates that nearest neighbors are considered only. The second term describes the screened Coulomb 
interaction. We denote with i?y = \R i — Rj\ the distance between lattice sites i and j and A -1 is the screening length 
which depends on the carrier concentration. The last term is due to the exchange and < ij > refers of pairs of ions 
in the same chain. While exchange plays a minor role in a Wigner crystal due to the large lattice constant it is of 
importance when crystallization takes place at high densities as it is the case here. 

Considering only the Coulomb term without taking screening into account except by using an effective charge, it 
has been demonstrated that the energy difference between a uniform distribution of holes and one with the holes 
concentrated in one family of chains is of the order of a few meV per formula unit [102]. A comparison with the 
effective / bandwidth of order 0.2 eV makes plausible that hole crystallization will take place. Note that there are also 
other charge-ordered configurations of low energy. For example, an ordering of charges in the form Yb 3+ — Yb 2+ — 




(4.1) 
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Yb 2+ — Yb 2+ — Yb 3+ ... in all four families of chains has also a low Coulomb- repulsion energy. In fact, von Schncring 
and Grin have found that the Coulomb-repulsion energy is slightly lower than the one with all Yb 3+ placed into one 
family of chains. However, that state is fourfold degenerate and therefore can further lower its energy by a Jahn- Teller 
distortion. This leads back to a state with short and long chains and therefore with holes concentrating in the short 
chains [30] . However, that order will not be a perfect one for the reason that the 4- fold degenerate state has a slightly 
lower repulsion energy. Thus, the optimal state must have incomplete order. Experimentally a 10 % deviation from 
perfect order is observed [22] . A different way of looking at imperfect charge order is to realize that the zero-point 
fluctuations of the holes on chains in the [111] direction necessarily leads to a spreading of holes into the long chains. 
For a rough estimate of the effect see Ref. [100]. 

The above physical picture justifies a description of the phase transition in terms of a band Jahn- Teller effect [30]. 
The corresponding Hamiltonian is 

4 4 

H=-tJ2 E «- a ^ + ^c-)-er^^A^ a +4 0) e2iV . (4.2) 

H=l <ij>a fJ.=l ia 

Here /j is a chain index and i labels a site on chain /i. Furthermore, er(r = T$) is a volume conserving strain 

order parameter which couples to the deformation potential A M = A[(<5 M i — (1 — <^i)/3] and is the associated 

elastic constant. N is the number of sites. A Jahn- Teller phase transition takes place if A 2 /(tep ^) > 3. From LDA 

calculations one can estimate At w 0.2eV. Furthermore, cp^/fi = 4 • 10 11 erg/cm 3 where D, is the volume of the 
unit cell. We choose A = 5eV which gives a Griineisen parameter of £l c = A/(4t) ~ 25 typical for intermediate 
valence compounds. With this set of parameters a transition temperature T$ ~ 250K is obtained in approximate 
agreement with the observed value [30]. Note that strong correlations not included in (4.2) may influence considerably 
the details of the transition [103]. As discussed in Sect. 5.3 the low-temperature specific heat of Yb^As^ is of the form 
C = -fT and therefore quite different from that of a Wigner crystal. This is due to the exchange interactions which 
are important in the crystalline phase. The charge ordering leads to quasi one-dimensional Heiscnbcrg spin chains 
and it is known that the latter have a specific heat linear in T. 

Yb 4 As 3 is most likely not the only material showing electron or hole crystallization. Other candidates are E114AS3 
and EU3S4 where Mossbauer and other measurements have shown that at low temperatures the Eu 3+ and Eu 2+ ions 
are at fixed lattice sites [130, 131]. More precisely, the valence fluctuation times must be larger than r = 10~ 8 sec, a 
typical testing time in a Mossbauer experiment. It should be mentioned that there has been an attempt to explain 
the experiments on E114AS3 and E113S4 by a Verwey transition [132 - 134]. 

V. HEAVY FERMIONS 

The investigation of metallic systems with heavy quasiparticle excitations has developed into an own branch of 
low-temperature physics. In most cases these systems contain Ce, Yb, U or Np ions as one of their constituents, 
implying that 4/ or 5/ electrons are essential. Examples are the metals CeAl 3 , CeCuiSii, CeRu^Sii, CeCu§, 
YbCu2Si2, UBe\3, UPt^, and NpBei^. For experimental reviews see Refs. [23,24] and [26,28], respectively. But also 
the electron-doped cuprate Nd,2- x Ce x Cu0 4 shows heavy-fermion behavior [32] in the range 0.1 < x < 0.2. Heavy 
quasiparticles have also been found in semimetals like ^64^53, Sm^Si or in some of the monopnictides and even in 
insulators like YbB\2 or SuiBq [104]. 

The following experimental findings define a heavy-fermion system: 

(a) A low temperature specific heat C — 7T with a 7 coefficient of order Umol^ 1 K^ 2 , rather than 
lmJmol~ 1 K~ 2 as, e.g., found for sodium metal; 

(b) A Pauli spin susceptibility xs which is similarly enhanced as 7; 

(c) A ratio R = n 2 k% xs/(^^lffl) (Sommerfeld- Wilson ratio) of order unity. 

Both quantities, 7 and xs are proportional to the quasiparticle density of states N*(0) per spin direction at the Fermi 
level. The latter is proportional to m*, the effective mass of the quasiparticles. When R is calculated the density of 
states 7V*(0) drops out. For free electrons R = 1, while in the presence of quasiparticle interactions R = (1 + Ffi) -1 
where Fq is a Landau parameter. When conditions (a)-(c) are met, one may assume a one-to-one correspondence 
between the low-energy excitations of the (complex) system like CeAl 3 and those of a free electron gas, provided a 
strongly renormalized effective mass m* is used and, in the case of semimetals or insulators, an effective charge e*, 
instead of the corresponding bare quantities. 
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Heavy-fermion behavior requires the presence of a characteristic low-energy scale in the system. The latter is 
usually characterized by a temperature T*. As the temperature T of the system exceeds T* the quasiparticles lose 
their heavy-mass character. The specific heat levels off, and the spin susceptibility changes from Pauli- to Curie-like 
behavior. With further increase of temperature the rare-earth or actinide ions behave more and more like ions with 
well-localized / electrons. 

A key problem is to understand the physical origin of the low-energy excitations. For a long time it was believed that 
the Kondo effect is the sole source of heavy quasiparticles. The physics of the Kondo effect is extensively described 
in a monograph [28] and a number of reviews [23-27]. However, by now it is known that also other effects may lead 
to heavy-fermion behavior. In all cases a lattice of 4/ or 5/ ions is involved though. In metallic systems this lattice 
couples to the conduction electrons. The latter are either weakly correlated like in CeAl 3 or strongly correlated like 
in the cuprates which may become high-T c superconductors. Strong correlations among the conduction electrons may 
influence substantially the physical properties of the system. Such a situation is encountered in Nd,2- x Ce x CuOi and 
as shown below it is here the Zeeman effect which is responsible for the formation of heavy-fermion excitations. In 
the semimetal Yb^As^ the heavy quasiparticles are intimately related to quasi one-dimensional chains of Yb 3+ ions 
which interact antifcrromagnetically with each other. It is well known that a Heisenberg chain has a linear specific 
heat C — jT at low temperatures and a Pauli like susceptibility. Thus, instead of having a single physical origin, 
heavy fermions may result from a variety of effects. 

Obviously, the low lying excitations, the main feature of heavy-fermion systems involve predominantly spin degrees 
of freedom. Evidential is the entropy S associated with the excess specific heat. It is of order S ~ ks^nvf per / 
site, where Vf is the degeneracy of the crystal-field ground state of the incomplete atomic / shell. In the following we 
discuss the three different routes to heavy-fermion behavior just outlined. It is likely that they will be supplemented 
by other ones in the future. 



V.l. Kondo lattices 



The essence of the single-site Kondo effect is the formation of a singlet ground state due to a weak hybridization 
of the incomplete 4/ shell with the conduction electrons. We derive the singlet wavefunction by starting from the 
Anderson impurity Hamiltonian 

H = J2<k)c+ n c km + e f J2n f m + j ]T n^n^, + (5.1) 

km m m^m' 

km 

Here /+ denotes the creation operator of an / electron in state m of the lowest J multiplet and — f^fm- The 
/-orbital energy is Cf and U is the / — / Coulomb repulsion. The c^ m create conduction electrons with momentum 
\k\ = k and the quantum numbers £ = 3, J and m. The hybridization between the / and conduction electrons is given 
by the matrix element V(k). Finally, H contains all those degrees of freedom of the conduction electrons which do 
not couple to the 4/ shell. The following ansatz for the singlet ground-state wave function is due to Varma and Yafet 
[105] 

|Vo>=A(l + ^Va(fc)/+c fem )|<A)> (5.2) 

where \4>o > represents the Fermi sea of the conduction electrons. It is closely related to the one suggested by 
Yoshida [106] for the ground state of the Kondo Hamiltonian. The variational parameters A and Aa(k) are obtained 
by minimizing the energy E =< V'ol-^lV'o > / < V'olV'o >• The latter is always lower than the one of the multiplet 
I ipm >= /ml^o >• The difference e is found to be 

e = -D exp[-\e f \l{ Vf N{Q)V 2 )] (5.3) 

and denotes the energy gain due to the formation of the singlet. Here D is half of the bandwidth of the conduction 
electrons and N(0) is their density of states per spin direction. It is customary to associate a temperature Tk, the 
Kondo temperature with this energy gain. The singlet-triplet excitation energy — e is often of the order of a few meV 
only, and provides a low-energy scale. When a lattice of / sites is considered instead of a single one, e.g., like in 
CeAl 3 the Anderson-lattice Hamiltonian is replacing Eq. (5.1). The energy scale ksTx is then replaced by a related 
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one, ksT* , which includes modifications due to the interactions between different / sites. The energy gain due to the 
formation of singlets competes with the one of magnetic /-sites interacting via the RKKY interaction [107]. The latter 
always wins for small enough hybridization V, because it is proportional to V whereas ksT* depends exponentially 
on V (sec Eq. (5.3)) and therefore is smaller in that limit. This seems to be the case in systems like CeAl 2 , CePbs 
and NpBeis which are antiferromagnets at low temperatures. 

In addition to T* there does exist another characteristic temperature T co h < T* below which the local singlet-triplet 
excitations lock together and form coherent quasiparticle excitations with large effective mass m* . The details of this 
transition are still an open problem, but de Haas-van Alphen measurements have demonstrated convincingly that at 
low temperatures the / electrons behave like delocalized electrons [108]. They contribute to the Fermi surface and to 
large effective mass anisotropies. It is surprising that one can calculate the Fermi surface of some of the heavy-fermion 
systems and determine the mass anisotropies with one adjustable parameter only. This is achieved by renormalized 
band-structure calculations (for reviews see [109, 110]). They are based on a description of the effective potential 
seen by a quasiparticle in terms of energy-dependent phase-shifts Tjf-{e) of the different atoms A. The index I refers 
here to the different angular momentum channels. As an example we discuss in the following the calculation of the 
Fermi surface of CeRuiSii [111]. The essential approximation is to use for the phase shifts the ones computed within 
the LDA, with the exception of the i = 3 phase shift of the Ce ions. Thus, only the r?£^ 3 (e) phase shift remains 
undetermined. It contains the strong correlations of the 4/ electrons and cannot be properly evaluated within the 
LDA. This approximation neglects virtual transitions between different crystal-field eigenstates caused by the coupling 
between conduction and 4/ electrons. (The mass enhancement of the conduction electrons in Pr metal falls into that 
category [112]). 

According to Hund's rules the ground state multiplet of Ce 3+ with a 4/ 1 configuration is j = 5/2. The multiplct 
j = 7/2 is sufficiently high in energy that it may be neglected and therefore J?^f 7 / 2 ( e F) = 0. Of the j = 5/2 multiplet, 
only the Kramers degenerate crystal-field ground state is taken into account, because it is the only one occupied at 
low temperatures. Therfore, near the Fermi energy only the phase shift function 7^? e (e)(T = 1,2) among the different 
I = 3 channels differs from zero. It contains the strong electron correlations and its form is unknown. In the spirit of 
Landau's Fermi-liquid theory we expand this function in the vicinity of e F and write 

ryf e (e) = if T e (e F ) + a(e - e F ) + 0((e - e F f). (5.4) 

The expansion contains the two unknown parameters r]^ e (e F ) and a. One of them, i.e., r]^ e (e F ) is fixed by the 
requirement that a Ce site contains one 4/ electron (n/ = 1). According to Friedel's sum rule this implies r^ e {e F ) = 
The remaining parameter a fixes the slope of the phase shift at e F . It therefore determines the density of states and 
with it the effective mass of the quasiparticles. We set a = (ksT*)^ 1 and determine T* by the requirement that 
the specific heat coefficient 7 calculated from the resulting quasiparticle dispersion agrees with the experimental 
one. Calculations of this form have explained and partially predicted [109,111] the Fermi surface and the large mass 
anisotropies of CeRu 2 Si2 (see Table 10 and also Fig. 10) [113, 114]. For more details on renormalized band theory 
we refer to the comprehensive reviews [109, 110]. 

When the temperature exceeds T co h the excitations lose their coherence properties and we are dealing with approx- 
imately independent scatterers. In that regime the specific heat contains large contributions from the incoherent part 
of the / electron excitations. 

The noncrossing approximation (NCA) is a valuable tool for treating the coupled 4/ and conduction electrons in 
that temperature regime [116-118]. It leads to a system of coupled equations of the form 

+00 

= £ / <kp m (<;)K + {z - <;) (5.5) 

m 

— OO 

E m (z) = - / -(*-?)■ 

J-00 

Here T = irN(0)V 2 and K±(z) are defined by 

—00 

where /(e) is the Fermi energy and N(e) is the energy-dependent conduction-electron density of states. The function 
T, a (z) and p a (z)(a — 0,m) relate to each other through 
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Pa(z) = ~-Im{R a (z)} (5.7) 

IT 

R a {z) = 



— e a Y, a (z) 

with e Q= o — 0, e a =m — £fm- The NCA equations have to be solved numerically [119]. However, one can find simple, 
approximate analytic solutions which have the advantage that crystal-field splittings can be explicitly included, a goal 
which has not been achieved yet by numerical methods. Once the p a (e) are known, one can determine, e.g., the 
temperature dependence of the /-electron occupancies n/ m =< /„/ m > through 

n fm (T) = j- J dep m (e)e-^-^ , (5.8) 

— OO 

where fi is the chemical potential and 

Zf= f p 1 e-^(R (z) + ^Rm(z)) (5.9) 
Jc 2m ^ 

is the partition function of the / electrons. Knowing the nf m (T) one can compute quantities like the temperature 
dependence of the quadrupolc moment of the / sites 

m 

The theory has been used to explain the observed Q(T) behavior of Xb in YbCu 2 Si2 [120, 121]. 

When T » T* , the / electrons can be treated as being localized. Via an exchange coupling they are weakly 
interacting with the conduction electrons and perturbation theory can be applied in order to study the resulting 
effects. 

A beautiful justification of the above scenario is the observed difference in the Fermi surface of CeRu2Si2 and 
CeRu-iGei which is shown in Fig. 10. When Si is replaced by Ge the distance between Ce and its nearest neighbors 
is increased. This causes a decrease in the hybridization of the 4/ electrons with the valence electrons of the neighboring 
sites. While in CeRuiSii the characteristic temperature is T* ~ 15.fr, it is practically zero in CeRu^Sii- De Haas-van 
Alphen experiments are performed at a temperature T ~ IK implying that for CeRu 2 Si2 it is T <C T* while for 
CeRu2Ge2 one is in the regime T >• T*. Therefore, the 4/ electron of Ce contributes to the volume enclosed by the 
Fermi surface of CeRu2Si2, but not of CeRu2Ge2- Indeed, Fig. 10 shows that the two Fermi surfaces have similar 
features, but the enclosed volumes differ by one electron. The Fermi surface of CeRu2Ge2 has a decreased electronic 
part and an increased hole part as compared with the one of CeRu2Si2- 



V.2. Zeeman scenario - Nd2-xCe x CuOA 



Low-temperature measurements of the specific heat and magnetic susceptibility have demonstrated the existence 
of heavy quasiparticlcs in the electron doped cuprate Nd2- x Ce x CuOi [32]. For x = 0.2 and temperatures T < 
IK the linear specific-heat coefficient is 7 = 4J/(mol ■ K 2 ). The magnetic susceptibility Xs is approximately T- 
indcpendent in that temperature regime and the Sommcrfcld- Wilson ratio is R ~ 1.8. The experimental findings are 
shown in Fig. 11. While these features agree with those of other heavy- fermion systems, there are also pronounced 
differences. In superconducting heavy-fermion systems like CeCu2Si2 or UPt^ the Cooper pairs are formed by the 
heavy quasiparticles. This is evidenced by the fact that the jump in the specific heat AC at the superconducting 
transition temperature T c is directly related to the large 7 coefficient, i.e., AC (T c ) / (■yT c ) w 2.4. The low-energy 
excitations are therefore strongly reduced below T c because one must overcome the binding energy of the pairs. But 
in superconducting Nd\.^CeQ.\^CuOi the formation of Cooper pairs has no noticeable effect on the heavy-fermion 
excitations. They remain uneffected by superconductivity 

A crucial difference between Nd2- x Ce x CuOi and, e.g. CeCw2<5'*2 are the strong electron correlations between 
the conduction electrons present in the former, but not in the latter material. In the two-dimensional Cu — O 
planes of N d2- x Ce x CuO^ with x > 0.1 we have to account for antifcrromagnctic fluctuations. There is considerable 
experimental evidence that these fluctuations are very slow at low temperatures. Consider undoped Nd^CuO^, an 
antiferromagnet with a Neel temperature of Tjy ~ 270 K. The exchange interactions between a Nd ion and its 
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nearest- neighbor Cu ions cancel because of the antiferromagnetic alignment of the Cu spins. Therefore one is left 
with the next-nearest neighbor Cu — Nd spin interaction. It is of the form as Cu S_ Nd and larger than the Nd — Nd 
interaction. The Schottky peak in the specific heat seen in Fig. 11 results from the spin flips of the Nd ions in the 
staggered effective field a < s Cu > set up by the Cu spins (Zccmann effect). It is also present in doped systems 
like Ndi.$Ceo.2Cu04 where antiferromagnetic long-range order is destroyed by doping. This can only be understood 
if the changes in the preferred direction of the Cu spins occur sufficiently slowly, i.e., slower than 10 -10 sec in the 
present case, to that the Nd spins can follow those changes adiabatically. Only then is a similar energy so that in 
Nd 2 CuOi required to flip a Nd spin. This physical picture has been confirmed by recent inelastic neutron-scattering 
and /iSR experiments [122, 123]. Spin-glass behaviour can be excluded. 

When Nd ions are replaced by Ce ions, the latter contribute approximately 0.5 electrons more to the Cu — O planes 
than the former. Thus a corresponding number of Cu sites are in a 3d 10 configuration. These sites have no spin 
and consequently they do not interact with the Nd ions. The extra electrons move freely in the Cu — O planes and 
therefore, the interaction of a Nd ion with the next-nearest Cu site is repeatedly turned off and on. It is this feature 
which results in heavy-quasiparticles. 

Two model descriptions have been advanced in order to explain the low-energy excitations of N di- x Ce x CuOi. 
One is based on a Hamiltonian in which the Nd — Cu interaction is treated by a hybridization between the Nd 4/ 
and Cu 3d orbitals. Usually it is much easier to extract heavy quasiparticles from such a Hamiltonian than from 
one with a spin-spin interaction like the Kondo Hamiltonian. The slow, antiferromagnetic fluctuations of the Cu 
spins are replaced by a static staggered field acting on them. This symmetry-breaking field also accounts for the 
strong correlations in the Cu — O planes because charge fluctuations between Cu sites are strongly reduced this way 
(unrestriced Hartree-Fock) . The Hamiltonian H reads therefore 

H = -t i a & a j* + h - c -) +hJ2 (Tei — ia t a i° ( 5 - n ) 

<ij>a icr 

+ V^a+fur + h.c.) + If ]T f+fi*. 

icr icr 

Here Q = (tt, tt) is a reciprocal lattice vector, denotes the positions of the Cu ions and h is the staggered field. 
The operators a+ , create an electron in the Cu 3d x i_ y i and Nd 4/ orbital, respectively. For simplicity, only one 
Nd site per Cu site is considered and one 4/ orbital with energy if is assumed instead of seven. The energies if and 
V are strongly renormalizcd quantities because of the 4/ electron correlations. 

The Hamiltonian (5.11) is easily diagonalized. Four bands are obtained, two of which are d-\ike (Cu) and two which 
are /-like (Nd). The dispersions of the four bands are given by 

E v {k) = ± lyJiekTif) 2 + W*, „=!,.. .,4 (5.12) 

where e& = (e§(fc) + /i 2 /4)^ and eo(fc) = —2t(cosk x + cosk y ). They are schematicly shown in Fig. 12. At half-filling 
only the lower / band is filled and the Schottky-peak contributions to C(T) arc due to transition from the filled 
lower to the empty upper / band. When the planes are doped with electrons the upper / band becomes partially 
filled resulting in low-energy intraband excitations with large effective mass. The latter follows from the quasiparticlc 
dispersion 

E »®* l '+&Tmy (5 ' 13) 

In the present case it is the Zeeman splitting of the / states which is responsible for the occurrence of heavy-electron 
behaviour. The effect of superconductivity on the heavy quasiparticles can be studied by adding an attractive part for 
the charge carriers in the Cu — O planes to the Hamiltonian. The latter can again be diagonalized. When the density 
of states is calculated one finds unchanged contributions from the / bands inside the BCS gap. They originate from 
the Nd spin degrees of freedom and explain why the heavy quasiparticles remain unaffected by superconductivity. 

The second model description of the Nd spins coupled to the Cu spin is based on stochastic forces acting on the 
latter [124]. They mimic the interaction of a Cu spin with the other Cu spins. In that case we start from the 
Hamiltonian 

H int =as Cu - S f , a>0 (5.14) 

describing the Nd — Cu interaction. For simplicity, both spins are assumed to be of magnitude S. We treat the 
vector Q = s Cu /S like a classical variable, subject to a stochastic force. We assume a Gauss-Markov process in which 
case the distribution function obeys a Fokker-Planck equation. The correlation function is then of the form 
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< £l(0)n(t) > = e~ Wrt (5.15) 

where D r can be obtained from the nonlinear a model [125]. Because there is no long-range order < Q(t) > = 0. 
The motion of the Nd spin is governed by the equation 

n(t) = Lu {n{t) x n(t)) (5.16) 
where n{t) = S_f/S and u>o = aS. The spectral function 

J(w) = [ dte lut < n(0)n(t) > (5.17) 

27T 7 

— oo 

is evaluated by making use of the corresponding stochastic Liouville equation. We find that I(u>) is of the form 

1 4D r 

T ^ = 3 W + (4£> r ) 2 + ( side P eaks at ^o). (5.18) 

While D r {T) vanishes as T —* in the presence of long-range order, it remains finite when the latter is destroyed 
by doping. A linear specific-heat contribution of the 4/ spin is obtained from 

oo 

d TT S(S + 1) f , 2 , r „, 

C(T W = - < ^ >= -L^ J d^ co JJ /2T) (5.19) 

o 

when D r (T = 0) 7^ 0. The side peaks of I(u>) give raise to a Schottky-type contribution. The calculated specific heat 
reproduces the experiments reasonably well. One shortcoming of the theory in its present form is the low-temperature 
spin susceptibility which follows from 

00 

Xim P (T) = ^(g» B ) 2 S(S + 1) / dJ-^tanh^. (5.20) 



Instead of a temperature independent Pauli paramagnetism we find Xim P (T) ~ £n(D r /T) at low T. This is possibly 
due to the neglect of Nd — Nd interactions. However, when evaluated for T = 0.4 K one obtains for Nd1.sCeo.2CuO 4 
a Sommcrfcld- Wilson ratio of R ~ 1.4. 



V.3. Hubbard route: F6 4 As 3 



The semimetal ^64^53 is an example of a system in which 4/ holes crystallize at low temperatures. The structure 
of the material and the transition from a metallic high-temperature phase to a semimetallic low-temperature phase 
were discussed in Sect. IV. Here we concentrate on the heavy-fermion properties at low temperatures. The following 
experimental observations are relevant in that respect. Measurements of the Hall constant demonstrate that at low 
temperatures only one carrier per 10 3 Yb ions remains. We interprete this being due to the Yb 3+ chains representing 
almost half-filled Hubbard systems. The resistivity is at low temperatures of the form p(T) = po + AT 2 like for a 
Fermi liquid. A linear specific heat is found at low T with a 7 coefficient of order 7 ~ 200mJ/(mol ■ K 2 ). The spin 
susceptibility is Pauli like and similarly enhanced as 7, giving rise to a Sommerf eld- Wilson ratio of order unity. No 
indication of magnetic order is found down to T = 0.045-ff, but below 2K the susceptibility starts to increase again 
indicating the presence of another low energy scale [126]. These findings strongly suggest heavy-fermion behavior 
which is further confirmed by the observation that the ratio A/j" with v = 2 is similar to that of other heavy-fermion 
systems. We reemphasize that the 7 coefficient exceeds the one, e.g., of Na metal by a factor of more than 10 2 despite 
a carrier concentration of only 1 per 10 3 Yb ions. This shows clearly that the low-energy scale must involve spin 
degrees of freedom of the Yb 3+ ions. Indeed, chains of antiferromagnetically coupled spins have a linear specific heat 
C = 7T. Inelastic neutron scattering experiments by Kohgi and coworkers [22] have demonstrated that the magnetic 
excitation spectrum is that of a Heisenberg chain with a coupling constant J = 25K. This spectrum leads to a 7 
value of the observed size. The physical origin of the heavy-fermion excitations is therefore very different here than 
that in the Kondo-lattice case. 
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A theory has been developed which explains consistently the above experimental findings [30]. It is based on 
interpreting the structural phase transition in terms of a band Jahn- Teller (CBJT) effect. This interpretation is 
suggested by the physical considerations outlined in Sec. IV in connection with the hole crystallization taking place. 
In the theory applied here the crystallization is due to a strong deformation-potential coupling which is quite common 
in mixed-valence systems. This potential has its origin in the Coulomb repulsion of the 4/ holes. 

The CBJT transition splits the fourfold degenerate quasi-ld density of states into a nondegenerate one correspond- 
ing to the short chains and a threefold one due to the long chains. The nondegenerate one is lower in energy and 
would be half filled if charge ordering were perfect and the holes were uncorrelated fermions. Instead, the holes are 
strongly correlated. Two holes on a site imply a 4/ 12 configuration for Yb and that has a much too high energy to 
occur. Therefore, we are dealing with an almost full lower (hole) Hubbard band rather than with an almost half-filled 
conduction band. Therefore, the system should be close to an insulator. That Yb i As 3 is a semimetal and not an 
insulator is most probably related to the nonvanishing hopping matrix elements between 4/ orbitals in the long and 
short chains. We have discussed in Sec. IV that the zero-point fluctuations of the 4/ holes lead to a partial (though 
small) transfer of holes to the long chains. We speak of self-doping if this transfer causes the gap in the excitation 
spectrum of (quasi) one-dimensional Hubbard chains to vanish. Accurate conditions for self-doping are not easily 
worked out, but a first step in this direction was done recently [127]. 

The phase transition is described by an effective Hamiltonian of the form of Eq. (4.2). 

With increasing charge ordering, correlations become more and more important because with the increase in con- 
centration of holes in the short chains their average distance decreases. Therefore, at low temperatures T the t — J 
Hamiltonian (see Eq. (3.2)) or a Hubbard Hamiltonian (3.1) must be used. Using the former and making use of a 
slave-boson mean-field approximation we arrive at an effective mass enhancement of the form 

to* _ t 
TO fc ~ tS + (3/4)x J 

Here nib denotes the band mass, \ = Xij =< faofjia > ; <5 is the deviation of the short chains denoted by 1 
from half filling and J = At 2 /U , where U is the on-site Coulomb repulsion between holes. With U = WeV one finds 
J = 1 • 1CT 3 eV and using x(T = 0) = (2/ir)sin(ir(l — 8)/2) with S = 10~ 3 one obtains a ratio of m* /m b ~ 100. 
The derivation of the mass enhancement hides somewhat the fact that spin degrees of freedom are responsible for the 
heavy quasiparticles. They become better visible when one sets S — 0, which is the case of no charge carriers. Even 
then one finds fermionic excitations with a large effective mass and corresponding heat coefficient 7. The theory was 
recently improved [103] by including in the Hamiltonian (5.21) an on-site Coulomb repulsion U between 4/ holes. 
This one-dimensional Jahn- Teller model can be solved exactly by a Lieb-Wu Bethe type ansatz. Of particular interest 
is that a self-doped distorted phase is obtained in a sizable regime of parameters. Since spin-wave-like excitations are 
responsible for the fermionic low-energy excitations associated with the specific heat and susceptibility at low T we 
are dealing here with charge-neutral heavy fermions in distinction to the charged heavy electrons, which appear, e.g., 
in CeAl 3 . Therefore, we speak of an uncharged or neutral heavy Fermi liquid. 

The physical interpretation given above allows for an explanation of another experiment. It has been previously 
found that an applied magnetic field of H — 4 Tesla has little influence on the 7 coefficient above 2 K, but suppresses 
7 considerably below 2 K [128]. This effect is unexpected, since one would have thought that the changes are of order 
(yU_B_ff/fcsT*) 2 and therefore very small. However, we can explain the experiments by providing for a weak coupling 
between parallel short chains. When linear spin-wave theory is applied, a ratio of order 10 -4 between interchain and 
intrachain coupling opens an anisotropy gap which modifies C(T) in accordance with observation [129]. 
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Figure Captions 



Fig. 1: Detailed presentation of one-body increments for NiO. From Ref. [61]. 
Fig. 2: Two-body oxygen-oxygen increments for MgO. From Ref [38]. 

Fig. 3: Van der Waals interaction in ionic crystals. Calculated values within the CCSD 
approximation versus values as obtained from London's equation (2.28). Circles: 
interactions between positive and negative ions; crosses: between negative ions. From 
Ref. [63]. 

Fig. 4: Contributions of different increments to the binding energy of MgO. From Ref [63]. 

Fig. 5: Schematic representation of the Cu Zd x i_ y i and O ^Vx{y) orbitals which are treated 
by a 3-band Hubbard Hamiltonian. 

Fig. 6: Spectral density of Ni as obtained from a 5-band Hubbard Hamiltonian (3.4 - 3.6) 
with U = 0.56, J = 0.22, A J = 0.031 (in units of the SCF bandwidth) when different 
approximations are applied: 

(a) full spectrum, 

(b) SCF approximation, 

(c) correlations included but with J = A J = 0, 

(d) when A J = 0, 

(e) assuming = 1. 
[From Ref [83]). 

Fig. 7: Spectral density of the Cu — O planes: 

(a) at half-filling, 

(b) for 25 % of hole doping. 

Oxygen and Cu contributions are shown by dashed and solid lines, respectively. Pa- 
rameter values are Ud — 8, U p = 3, t pp = 0.5, e p — e<j = 4 in units of t p d- Note, 
that with hole doping spectral weight is shifted from the upper Hubbard band to the 
region close to the Fermi energy (dotted line). (From Ref. [87]). 

Fig. 8: Spectral density obtained by diagonalization of a cluster of (Cu02)i'- 

(a) at half-filling, 

(b) for 25 % of hole doping. 

Parameter values are similar to those in Fig. 7, i.e., Ud — 8.8, U p — 4.24, t pp — 0.41, 
£d = 3.37. Note the similarities with Fig. 7. S denotes the Zhang- Rice singlet and 
T indicates the spin-triplet contribution. Data has been broadened by a lincwidth. 
From Ref. [90]. 

Fig. 9: (a) Structure of Yb 4 As 3 , 

(b) four families of chains on which the Yb ions are located. 

Fig. 10: Fermi surfaces for CeRu^Sii and CeRu2Ge2 as derived from de Haas - van Alphcn 
measurements [113 - 115]. The upper part is due to holes and the lower part is due to 
electrons. In CeRuiGei the 4/ electron of Ce is well localized while in CeRuiSi2 is 
participates in the formation of the Fermi surface. The volume enclosed by the Fermi 
surface differs therefore by one electron resulting in an increase of the hole part and 
a decrease of the electron part in CeRu2Ge2- 

Fig. 11: Observation of heavy-fermion excitations in Nd2- x Ce x CuOi. 

(a) specific heat C P (T), 

(b) C P (T)/T, 

(c) Spin susceptibility for an overdoped sample with x = 0.2. 
From Ref. [32]. 

Fig. 12: Schematic plot of the quasiparticle bands of N d,2- x C e x CuO ^ for i ^ 0. The Fermi 
energy is indicated by a dotted line. Solid lines: /-like excitations, and dashed lines: 
d-like excitations. 
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Table Captions 



Table 1: 



Table 2: 
Table 3: 

Table 4: 

Table 5: 

Table 6: 

Table 7: 
Table 8: 



Cohesive energy per unit cell (in eV) of the elemental semiconductors in SCF approxi- 
mation (Eg°Q F ) and with inclusion of correlations (E^ r ). Percentages in parenthesis 
are with reference to the experimental values (E^p) which include an estimate of 
the atomic zero-point fluctuations [47]. For comparison results of other methods are 
shown: local density approximation (E^ A ) [48, 49]; generalized gradient approxima- 
tion (Eg£ A ) [50]; quantum Monte Carlo (£^ c ) [51]. From Rcf. [37] with corrections 
included due to atomic spin-orbit splittings. 

Cohesive energy (in eV) of the III - V semiconducting compounds. The notations are 
the same as in Table 1. From Ref. [37]. 

Lattice constants for the group IV semiconductors in Angstrom: in SCF approxima- 
tion (oscf); including a core polarization potential (a cpp ); with inclusion of corre- 
lations (a corr ); Percentages in parenthesis are with reference to experimental values 
(a exp ). One notices a strong influence of core polarization. From Ref. [37]. 

Bulk modulus of the group IV semiconductors in Mbar in SCF approximation (Bscf)', 
including a core polarization potential (B cpp ); with correlations included (B corr ); 
Percentages in parenthesis are with reference to experimental values (B exp ). From 
Rcf. [37]. 

Parameters in the BOA as obtained by fitting the results for a single and for neigh- 
boring bonds to those of the ab initio calculation with a minimal basis set. 8E % J^F 
denotes those contributions which come from more distant than nearest-neighbor 
bonds. From Ref. [54]. 

Cohesive energy (in eV) of three oxides. The notation is the same as in Table 1. E C °J^ 
from Ref. [56], Efg F from Ref. [57] (MgO), [58] (CaO) and [59] (NiO). Calculations 
for -E^orr are 0n a CCSD(T) level (MgO, CaO) and by applying quasidegenerate 
variational perturbation theory (QDVPT) [60]. 

Lattice constants (in A) of three oxides. Same notation as in Table 3. For ascF see 
Refs. [57 - 59]. 

One-body increments for MgO, CaO and NiO in eV. From Refs. [38, 40, 61]. 



Table 9: Sums of different two-body increments in eV. The last line gives the sum of all 2-body 
increments. From Refs. [38, 40, 61]. 

Table 10: Measured [114] and calculated [111] mass ratios for CeRu2Si2- Shown are extremal 
areas of the Fermi surface (in megagauss) and the effective mass ratios m* /mo- Unlike 
the LDA, the rcnormalized band theory (RB) reproduces the large measured mass 
anisotropies. 
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Si 


Ge 


a-Sn 


EfjCF 
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6.18 


4.25 


3.65 




(71%) 


(66%) 


(53%) 


(53%) 


E corr 


14.36 


8.84 


7.02 


6.13 




(95%) 


(94%) 


(88%) 


(90%) 


rpcoh 
-^cxp 


15.10 


9.39 


8.00 


6.83 


JTTCOll 

fi LDA 
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10.59 


9.06 




rncoh 
i'GGA 




8.79 


6.83 




TTICOh 

"QMC 


14.90 


9.76 
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BN 


BP 


BAs 


A1P 


AlAs 


AlSb 


TTlCOh 

-C'SCF 
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5.39 


4.71 


3.97 




(67%) 


(60%) 




(64%) 


(59%) 


(60%) 


y^rCOll 
_C/ corr 


12.38 


9.36 


8.57 


7.95 


7.18 


6.31 




(91%) 


(90%) 




(95%) 


(90%) 


(95%) 


rpCOh 
-^cxp 


13.61 


10.39 
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8.00 


6.61 




GaP 


GaAs 


GaSb 


InP 


InAs 


InSb 


rpcoh 


4.00 


3.54 


2.97 


3.86 


3.51 


3.18 




(53%) 


(53%) 


(49%) 


(57%) 


(54%) 


(55%) 


_c/ corr 


6.69 


6.20 


5.39 


6.37 


5.96 


5.47 




(91%) 


(93%) 


(89%) 


(94%) 


(92%) 


(94%) 


£££ 


7.37 


6.69 


6.07 


6.80 


6.50 


5.80 
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Si 


Ge 


Sn 


ascF 


3.5590 


5.4993 


5.7516 


6.6001 




(-0.2%) 


(+1.2%) 


(+1.7%) 


(+1.7%) 


Q-cpp 




5.4662 


5.6653 


6.4549 






(+0.6%) 


(+0.2%) 


(-0.5%) 


Clcorr 


3.5833 


5.4256 


5.6413 


6.4443 




(+0.5%) 


(-0.1%) 


(-0.3%) 


(-0.7%) 




3.5657 


5.4317 


5.6575 


6.4892 
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Si 


Ge 


Sn 


#SCF 
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1.038 


0.961 


0.638 




(+9%) 


(+5%) 


(+31%) 


(+20%) 


B C pp 




1.009 


0.889 


0.562 






(+2%) 


(+21%) 


(+6%) 
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0.979 


0.711 


0.510 




(-5%) 


(-1%) 
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(-4%) 


^cxp 


4.42 


0.99 
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0.531 
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Sn 
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13.8 
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3.20 


V 
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1.44 


1.25 


Vi 
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0.45 


0.38 


EXf (BOA) 
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-1.19 


-0.91 


-0.83 


S E mtor 


-0.45 


-0.42 


-0.40 


-0.38 
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MgO 


CaO 


NiO 


&SCF 


5.85 


5.92 


4.37 


rpcoll 


7.82 


8.06 


7.00 




(96 %) 


(93 %) 


(93 %) 


rpcoh. 
^exp 


8.14 


8.65 


7.49 
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MgO 


CaO 


NiO 


ascF 


4.191 


4.864 


4.264 


Q>corr 


4.184 


4.801 


4.164 


&exp 


4.207 


4.803 


4.170 
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MgO 


CaO 


NiO 


lattice const. 


4.21 


4.81 


4.17 


o -> o A ~ 


-2.04 


-2.05 


-2.14 


X -> X' z+ 


0.99 


1.00 


1.58 




-1.05 


-1.05 


-0.56 



Table 8 
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MgO 


CaO 


NiO 


o-o 


-0.36 


-0.15 


-0.56 


X-O 


-0.41 


-0.78 


-1.56 


X-X 


-0.002 


-0.02 


-0.04 




-0.77 


-0.96 


-2.15 
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